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Abstract 

It was recently proposed by the second author to consider lattice formulations of QCD 
in which complete actions, including the gauge part, are built explicitly from a given 
Dirac operator D. In a simple example of such theory, the gauge action is proportional 
to the trace of Ginsparg- Wilson operator D chosen to define the quark dynamics. This 
construction relies on the proposition that the classical limit of lattice gauge operator 
tv D(x,x) is proportional to tiF 2 (x) (up to an additive constant). Here we show this 
Q ■ for the case of the overlap Dirac operator using both analytical and numerical methods. 

tIh ■ We carry out the same analysis also for the tensor component of D, which is similarly 

related to the field-strength tensor F, and obtain results identical to our previous 
CN ■ derivation that used different approach. The corresponding proportionality constants 

are computed to high precision for wide range of the negative mass parameter values, 
■<^j- | and it is verified that they are the same in finite and infinite volumes. 

t> 

(N 

o 

oo : 1 Introduction 
o 

To define various versions of lattice QCD (LQCD), one usually follows a standard route 
^ ■ where the gauge part of the total continuum action is treated independently from the part 
involving fermions. Indeed, one normally writes down a candidate for lattice gauge action as 
an explicit function of link variables with appropriate symmetries, such that its continuum 
limit on smooth backgrounds coincides with the continuum expression. On the other hand, 
the basic entity in the construction of fermionic part is the lattice Dirac operator D, itself 
a collection of functions in gauge variables (matrix elements). In a standard treatment, one 
does not attempt to functionally relate gauge action to matrix elements of D. 

It has recently been suggested by one of us [H [2] that exploring lattice gauge theories 
with gauge and fermionic parts of the action explicitly related (coherent LQCD) could be 
beneficial for studies of QCD vacuum structure. This is particularly attractive if D is a 
chirally symmetric operator of Ginsparg- Wilson type. Q One unusual feature of coherent 



■"■It is worth emphasizing that the notion of coherent LQCD can be well-defined even if D is not chirally 
symmetric, assuming that a suitable function f(D) is used to define required gauge operators. For example, 
the choice f(D) = D A produces coherent LQCD for arbitrary lattice Dirac operator D Q]. 
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formulations is the possibility of incorporating explicit relations between gauge and fermionic 
aspects of the theory. For example, it was suggested that QCD with Nf quark flavors can 
be regularized via lattice formulation in which the usual gluon kinetic term is traded for 
additional quark flavor(s) whose mass controlls the gauge coupling (symmetric logarithmic 
LQCD) [HE]- It remains to be seen if one can gain theoretical or computational advantage 
by introducing interrelations of this type into the lattice-regularized theory. 

The validity of various coherent versions of LQCD hinges on conjectures of locality and 
proper classical limit for scalar and pseudoscalar gauge densities associated with operator 
functions f(D) used in the construction. If D is a Ginsparg-Wilson operator, which we 
assume from now on, then these are expected to be technically non-trivial issues due to 
the fact that D is necessarily a non-ultralocal operator in fermionic variables [3], and its 
dependence on gauge degrees of freedom is also expected to be generically non-ultralocal. 

In this paper, we discuss the issue of classical limit for coherent LQCD, where the gauge 
action is based on Tr D. More precisely, we will focus on standard overlap operators [I] 
constructed from Wilson-Dirac matrix (negative mass — p, Wilson parameter r), and show 
that the classical limit of tr D { - P ' r \x ) x) is proportional to tr F /lu (x)F fliy (x) up to an additive 
constant [T]. It should be emphasized that in formulations with f(D) = D one does not 
have to deal with the issue of locality. Indeed, the required locality properties follow from 
locality of D, which is imposed to begin with. For overlap Dirac operators the locality was 
studied e.g. in Refs. [5j E], and the specific issue of locality in gauge variables was indirectly 
checked via reflection positivity considerations in Ref. [7j. 

In addition to justifying simple coherent LQCD as a valid regulator, the validity of con- 
jectured classical limit for tr D^ p ' r \x, x) will also make possible its use as a coherent scalar 
partner to overlap-based pseudoscalar density (topological density) [H [H [101 [11] . The asso- 
ciated noise reduction due to non-ultralocality is expected to be useful for obtaining refined 
information on the fundamental vacuum structure first discovered in overlap-based topolog- 
ical density [121 [13], as we U as for conventional uses such as studies of glueball spectrum [13] . 
Needless to say, for these purposes, it would in fact be desirable to have all interesting com- 
posite fields defined coherently. Important example is the field strength tensor F^, which 
can be constructed from tensor component of D( p,r >{x,x) [TBI, [Tj IT5]. Direct derivation has 
shown that in this case the required classical limit is indeed correct [TH1 [T7] . 

In principle, one could employ methods of Ref. [17] to derive the classical limit also in 
the scalar case considered here. However, the corresponding calculation is technically rather 
involved, and we thus take a different approach. In particular, the classical limit is first calcu- 
lated for the class of constant gauge fields in infinite volume. This derivation is rigorous (even 
though technically it relies on a numerical evaluation of certain Brillouin zone integrals), and 
allows us to determine the corresponding proportionality constants c s (p,r) in a comparably 
simple calculation. Next, we show that our arguments also apply for constant fields in finite 
volume, at least for standard periodic and antiperiodic boundary conditions. Finally, we 
perform a direct numerical evaluation of the classical limit on generic non-constant back- 
grounds and find identical results again. It is thus safe to conclude that the conjectured 
classical limit is indeed realized. Simultaneously with the scalar case, we apply our methods 
also to the tensor spinorial component and reproduce the results reported in Ref. [T7J, thus 
cross-checking both approaches. The proportionality constants c s , c T and scalar subtraction 
constants relevant for practical applications are given. 
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2 Formulation of the Problem 



Before we formulate the statements that we will be concerned with in this paper, we wish 
to fix our notation and specify few conventions. 

(i) Continuum Gauge Fields. The SU(3) gauge field A^x) (traceless Hermitian 3x3 
matrices) is related to the field-strength tensor via d 

F^(x) = dpA v (x) -d v A^x) +i[A lx (x),A u (x)} (1) 

while the covariant derivative acts as 

D^{x) = (d^ + iA^x))^) [D^D^x) = iF^x) <f>(x) (2) 

(ii) Classical Fields. By classical continuum fields on R 4 (or on a torus) we mean gauge po- 
tentials A^ix) smooth (differentiable arbitrarily many times) almost everywhere. If classical 
backgrounds contain singularities, the classical continuum limits are assumed to be taken at 
its non-singular points. While our conclusion is expected to be valid for all classical fields 
(due to the locality of operators involved), in what follows we will only consider classical 
fields that are smooth everywhere. This avoids various inessential technical complications 
related to transcribing singular fields on to the lattice. 

(Hi) Transcription of Classical Fields to Hypercubic Lattice. Hypercubic lattice is superim- 
posed on R 4 via the correspondence x = na (x £ R 4 , n £ Z 4 ), where a is the lattice spacing. 
Smooth continuum field is then transcribed to the lattice field via 



U n ^{a) = Vexp^ia J ds A^an + (1 — s)afi)J 



(3) 



where V is the path ordering symbol and fi is a unit vector in direction \i. 
(iv) Overlap Operators. Standard overlap Dirac operators are defined by [I] 

D M =p[l + X(X ] X)-?] X = D w -p p£(0,2r) r>0 (4) 

where Dw = 4rl — \K is the massless Wilson-Dirac operator and K, the Wilson's hopping 
matrix with Wilson parameter r, is given by 

■Kn,m ^ (j" T/x) 0~ri+n,m ~\~ {j* T/-t) ^n— )J,,fJ, ^n—[i,m (5) 

For calculations in this paper, it is convenient to work with the rescaled matrix X namely 

X — ► X = 2kX = I — kK k= (6) 

8r-2p y 1 

The form of the overlap operator (jlj) in terms of X is identical to that in terms of X. 



2 Note that this differs from conventions of Ref. pQ, where anti-Hermitian gauge potentials were used 
instead. The equations used here can be obtained from equations of Ref. pQ via substitutions An(x) — > 
iA^ix), F^ v (x) —> iFfj, v (x). The values of constants c s , c T are the same in both conventions. 
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(v) Classical Limit Conjecture. Our main goal in this paper is to support the Conjecture C3 
of Ref. [1] for the specific case of standard overlap Dirac operators. In the infinite volume, 
one can state this as follows. H 

Conjecture 1. Let A^(x) be arbitrary smooth SU(3) gauge potentials on M 4 . If U(a) = 
{ U n ^(a) } is the transcription of this field to the hypercubic lattice with classical lattice 
spacing a, and I = { U n ^ — > I c } is the free gauge configuration then 

tv sc ^D ofi (U(a)) - A,,oCQ) = c 5 a 4 tr c F^(0)F^(0) + 0(a 6 ) (7) 

where D = _D^> r ) is the overlap Dirac operator and F^ v (x) is the field-strength tensor asso- 
ciated with Afj,(x). The constant c s = c s (p,r) is non-zero and independent of A^(x) at fixed 
p and r. 

We will also discuss classical limits for tensor components of D^ p ' r >. The corresponding 
statement is obtained from the above by replacing Eq. (J7J) with pp, [17] 

tr a <vA,, (l7(fl)) = c T a 2 F^(0) + 0(a 4 ) (8) 

where = ^[TahTJ-u 

We emphasize that the above statements (with identical proportionality constants c s (p,r), 
c T (p,r)) are also expected to be true in a finite volume. We will verify this here for the 
practically relevant case of a torus. In what follows, we will frequently skip explicitly de- 
noting dependences of operators/constants on parameters p and r, but they are implicitly 
understood. 

3 Constant Fields in Infinite Volume 

In this section we establish that the leading term of Eq. (j7J is indeed correct for the class of 
constant gauge fields A^(x) = A^. This is the simplest class of classical fields leading to non- 
zero field-strength tensor (for non-Abelian gauge fields), while the resulting simplifications 
make the calculation reasonably manageable. 

For constant classical backgrounds the transcribed lattice field at lattice spacing a has 
the simple form 

U n ,n( a ) = exp(mA M ) (9) 

and the overlap Dirac matrix is translation invariant D n ^ m = D n ^ m ^. The calculation of 
classical limit simplifies in the Fourier space, where the operator is diagonal in space-time 

3 The subscripts "s", "c" in expressions that follow refer to "spin" and "color" respectively. Thus, for 
example, I c denotes the identity matrix in color space (3x3 matrix) and tr sc denotes the trace in both spin 
and color. 

4 We use Hermitian 7-matrices 7^ = 7^ throughout this paper. 

It is worth mentioning at this point that one can use perturbation theory techniques to obtain expansions 
needed [TUl HS1 HH1 HO] • Here we proceed without reference to weak coupling perturbation theory, as done 
in [TTJCE7]. 
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indices. For arbitrary translation invariant operator 0„ im we define the diagonal Fourier 
image O(k) via 

(k) = Y^e- %{n - m)k O n ^ m O n , m = ^4 J d^ke^-^ k O(k) (10) 

where the integration over momentum variables runs through the Brillouin zone. The defini- 
tion of 0(k) requires the convergence of the above infinite sum, for which it is sufficient that 
O is local. For arbitrary constant gauge field the overlap Dirac operator D in transcribed 
lattice background U^a) is guaranteed to be local [5] at sufficiently small classical lattice 
spacing a, and thus its Fourier transform is well defined. 

To evaluate D(k), we start from the expression for D in terms of matrix X (see Eq. (J6])) 

-D = 1 + X—^L= (11) 

p Vjox 

The Fourier transform X(k) of X can be found straightforwardly and is given by 

X(k) = Igc — 2kt cos( aA^ + k^) + i2/t^^7 M sin^aA M + k^j (12) 

Moreover, the Hermitian matrix B = X^X is strictly positive-definite for sufficiently small 
a in arbitrary constant background [5]. Consequently, the Fourier image of B" 1 ^ 2 is simply 
given by (X^k)X(k))- 1 / 2 , and thus 

-D(k) = l sc + X(k)— r 1 (13) 

p K ' K Wx^k)x(k) y } 

3.1 Expansion in Lattice Spacing 

For the purposes of Conjecture 1, we are interested in the Taylor expansion (in a) of 

-D 00 = -—^-r / d 4 kD(k) = I sc + — ^—r [ d 4 kX(k)^ 1 (14) 
P 0,0 P (2vr) 4 J K } (2tt) 4 J K } ^X^k)X(k) K ' 

which can be obtained term by term from the Taylor expansion of X(fc)5~ 1/2 (£0 due to its 
analyticity in both k and a for sufficiently small a. 

To fix the notation that we will use for various operators, we generically write 

oo 

0{k,a) = ^O n {k)a n = O {k) + 5O{k,a) (15) 

n=0 

where the dependence of O(k) on a will usually not be shown explicitly. The leading terms 
for X(k) and B{k) are given by0 

cosfc M + i 2k yj 7 M sin k^ (16) 

H (i 



6 The gap in the spectrum of B tends to zero when excluded boundary values p = and p = 2r are 
approached. 

7 Note that in the formulas that follow we denote I sc simply as "1" with spin-color structure understood. 
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and 

B (k) = (l-2Kr^2 cos k^j + 4/t 2 ^ sin 2 A; M (17) 

A 1 A* 

There are two points to note here. (^j B (k) is proportional to identity and is positive- 
definite for all k, except when nr = 1/8, 1/4 which corresponds to excluded boundary values 
of the p-parameter 0, 2r. This means that B (k) can be commuted and inverted freely 
in the expressions that follow, (ii) The matrices SX(k) and SB[k) have norms that can 
be made arbitrarily small by lowering the lattice spacing sufficiently, and can thus serve as 
"perturbations" in matrix expansions around X (k) and B (k). 

The expansion of X(k)B- 1 ' 2 (k) in lattice spacing, needed in (iT4"|) . can now proceed via 
expanding X(k) and B~ l l 2 {k) separately and then combining the results. In case of X(k), 
the expansion is obtained directly from (fl2l) . while for B^ 1 ^ 2 (k) we use 

= (B + 8B)-? = B~ h (l + 8B)~? SB = B$ l SB (18) 

Properties (%) and (ii) guarantee that 5B is a valid perturbation for the expansion of the 
inverse square root, yielding 



B~? = B 7 



Note that we have to keep all the terms up to the 4-th order since we are eventually interested 
in 0(a 4 ) contribution and SB is 0(a). Indeed, the explicit formula for B(k) in terms of its 
Clifford decomposition is given by 

B(k) = I s x $(k) + 7a , x V M (fc) + x 7^{k) (20) 

where 

§(k) = (l-2KrJ2 C0S ( aA n + k n)) + ^ sin ( a \ + k n? ( 21 ) 

A* A* 

V M (A;) = iAk 2 [sin(aA M + cos(aA, + fc„)] (22) 

V 

7^ v (k) = i2K 2 ^sin^aA^ + k fl ),sin(aA u + k u ) (23) 

and [, ] denotes the commutator. From this one can inspect directly that while V M and 7^ v 
are 0(a 2 ), the scalar part contains both the constant and the linear term, and thus SB is 
indeed 0(a). 

The equations of this subsection together with Eq. fTl2|) define the procedure of evaluating 
various classical limits completely. The calculation is straightforward, but turns out to be 
technically still quite involved if one wants to obtain a complete expansion of D(k) (and thus 
of A),o) all the wa Y U P to order a 4 . 
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3.2 Sample Computation — First Order 

To see the characteristic nature of the computations involved, let us now expand D(k) up to 
order a. For this purpose, we will need expansions of X(k) and B(k) up to 0(a). We find 

X(k) = X (k) + a V 2n(rsmk fl + i^cosk^) + 0(a 2 ) (24) 

, V » ' 

and 

B(k) = Bo(k) + a'^^A^ l 4«r(l — 2nr cos k v ) sin k^ + 8k 2 sin k^ cos k^ +0(a 2 ) 



(25) 

Note that f^{k) comes entirely from S(k). Using the expansion ffl9l) up to 0(5B) we then 
obtain D(k) = D (k) + aD x {k) + C(a 2 ) where 

- D (k) = 1 + 5 (A;)- 1/2 X (A;) (26) 



- Di(k) = J2 A » B °(k)- 3/2 \g»(k)B (k) - \f^{k)X {k)] (27) 



and 

P 

The leading term D (k) is the Fourier transform of the free overlap operator as expected. 

Let us now focus on the first-order term which will determine (after inverse Fourier 
transform) the first-order term of D 00 . Using transformation properties, one can argue 
that this first-order term should in fact vanish. Indeed, Dq q is a gauge covariant operator 
and this covariance has to be preserved order by order in the Taylor expansion. Thus, the 
coefficient of the linear term is expected to be a dimension one (continuum) gauge covariant 
operator (spin-color function of the gauge field), but such operator does not exist. On the 
other hand, the result (1271) shows that Di(k) is not identically zero. This does not imply 
any contradiction as long as the associated inverse Fourier expression vanishes. One can see 
that this is indeed the case by evaluating /i M (/c) of Eq. fl27]) explicitly, yielding 

h,{k) = rf{k) + vyphfo) (28) 

where h^(k) is antisymmetric with respect to k^ — > — k^, while h s ^{k) is symmetric, fully 
diagonal and can be written SIS db derivative, namely @ 

h%k) = 2k A sin k^B'^k) (29) 
Consequently, J d A k Di(k) = and we have 

A,, (tf(a)) = A),o(I) + 0{a 2 ) (30) 

While relatively simple, this calculation fully illustrates the issues encountered also at higher 
orders, where the number of terms to deal with grows rather quickly. 



3 We thank the referee for pointing this out to us. 
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3.3 Tensor Component and the Second Order 

The tensor part of _Do,o first appears at second order in classical lattice spacing. Using 
the technique described above we straightforwardly arrive (see also Appendix [A]) at the 
expression in Fourier space, namely 

>> = £[i^^K*^ + -" < 31 > 

where F^ u = i [A^, A u ], a^ v = h [7^, 7„], and ". . ." represents terms that are either odd in k^ 
or are partial derivatives of analytic functions with respect to thus not contributing upon 
transition back from Fourier space. The expression in the bracket is a real-valued function 
of the momentum and it should be pointed out that from now on we use the symbol B (k) to 
interchangeably denote a color-spin matrix proportional to identity, as well as its diagonal 
matrix element. The function t^k) is given by 

tfiu{k) = sin 2 k^ cos k v + sin 2 k v cos k^ — ( 4 — — — cos k a J cos fc M cos k v (32) 

a 

The constant c T = c T (p,r) of Eq. (|Sj) is then specified by 

c T( p r) = 2r P J_ [ d 4 k-^- (33) 

and doesn't depend on p,u due to hypercubic symmetries. One can easily check that this 
result is identical to the one we obtained previously in Ref. [T7] . 

The parameter dependence of c T is more transparent if one introduces a rescaled mass 
parameter 

P = ~ P G (0, 2) (34) 

r 

and the rescaled B (k), namely 

*(*) = J^2 B ^ k ) = (4r-p) 2 B (k) (35) 



With that we finally have 



^ = ^ 1 (S Ml < 3e > 



where 

= ^ sin 2 (Ag + r 2 ( 4 - p - J] cos A; M ) 2 (37) 

A 4 A* 

depends on both p and r, while 

tfiu(k) = sin 2 /c^ cos ac„ + sin 2 fc„ cos A; M — (4 — p — cos k a ) cos k^ cos A;^ (38) 




I , , , , , , , u 1e-12 1 1 1 1 1 1 1 1 ■- 

5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 

N N 



Figure 1: Riemann sums for proportionality constant c T of the tensor term. These results 
were computed at r = 1 and k — 0.19 (p = 26/19). 



only depends on p. 

For explicit evaluation of c T we study the convergence of Riemann sums associated with 
the corresponding Brillouin zone integral. More precisely, we partition the integration do- 
main [0, 27r] 4 into N 4 cubes of volume (27r/iV) 4 = (5k) A and define the corresponding Riemann 
sum c T [N] via 

1 ^TfAri ^ ST^fxu\^ tp V {k) 1 tp V {h) 



2 p r 2- " ( 2vr )4 2^°^ z 3/2(jfe) " 7V4Z^ 2 3/2 (A .) 



A- 



where the discrete momenta associated with elementary cubes have components 

k =/W" for/i=l,2,3 
^ 1(^)^ + 1/2) for/i = 4 U' 1 '---^ l 4U ^ 

Note that this choice of discrete momenta exactly corresponds to the situation on the latti- 
cized torus of size N with mixed periodic (spatial) and antiperiodic (time) boundary condi- 
tions, which is frequently a preferred setup in finite volume. □ Due to the asymmetric choice 
of the discretization, the value of c T [iV] at finite iV will depend on whether the pair /i, v 
chosen for evaluation is of space-space or space-time kind, but the difference has to vanish 
in the N —>■ oo limit. In Fig. [T] we show the dependence of c T [iV] on iV for both cases, 
showing that they converge exponentially in iV to the common non-zero limit. The value 
labeled as "exact" includes only digits determined to be stabilized under increasing N in the 
calculation. 



3.4 Scalar Component and the Fourth Order 

Our main goal in this paper is the evaluation of classical limit for the scalar spinorial compo- 
nent of D 0)0 , and its complete trace in particular. Since the expectation is that the leading 
term (up to a constant to be subtracted) will appear at the 4-th order, one has to carry out 

9 The rationale for this is that it will allow us to argue that, for constant fields, our results on infinite 
lattice can be mapped exactly to this specific finite-volume case. 
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expansions described in Sees. 13. 1[ 13.21 up to that order. Algebraic manipulations involved in 
this are rather extensive and we have used Mathematica for required symbolic manipulations. 
The result can be written in the form 

tr sc D , {U{a)) = tr sc D , {I) + a 4 J 4 + 0(a 6 ) (41) 

The odd powers are absent in the scalar spinorial component (do not contribute upon taking 
the spinorial trace of Dq$). The quartic term is given by 

h = ^U^ticAl + ^2 ti c A„A„A V A V + ^ (3^ tr c A^A„A^A„ (42) 

To describe the results of calculation for constants i^, a^, j3^ u we now have to intro- 
duce few conventions. All the coefficients to compute will be defined by functions in the 
momentum space of the form 

J, c) = -i— V c(n) cos(hr cos(A; 2 )" 2 cos^)" 3 cos^)" 4 (43) 

where n = (ni, n-2, n^, n 4 ) is a four-component vector of non-negative integers, j is a single 
half-integer and c(n) is real-valued. Each function of interest is thus specified by j and the 
finite list of non-zero coefficients c(n). 

To reduce the amount of information needed for evaluation of the coefficients, we will use 
the fact that, by means of the inverse Fourier transform, it is only the mean value of F that 
is relevant. More precisely, we are interested only in 



(44) 



where we introduced the notation I(j,n) for integrals of the form 



f d 4 k cos(fc 1 ) Wl cos(fc 2 ) n2 cos(fc 3 ) n3 cos(fc 4 ) n4 

Since n) is completely symmetric with respect to indices n M , we can restrict the sum over 
n in (1441) only to n such that n\ > n 2 > n 3 > n 4 , and redefine the coefficients c(n) — ■> c(n) 
accordingly. We then have 

F{j,c)= C H 7 0» ( 46 ) 

ni>rt2>n^>n4 

and we specify our results by listing j and c(n) for each case in question. 

The last piece of convention we need to specify concerns the choice of discretization for 
performing Riemann sums in the actual evaluation of momentum integrals. While final re- 
sults have to be independent of that choice, our discussion of classical limits in the finite 
volume (see Sec. H]) can be mapped directly on the discussion here if Riemann sums are de- 
fined in a particular manner. Moreover, consistency over different discretizations represents 
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Figure 2: Riemann sums for coefficients of Eq. ( )42l) . These results were computed at 
r = 1 and k = 0.19 (p = 26/19). 



an additional check on our results. In what follows, we will refer to the discretization defined 
by Eq. (j40p as antiperiodic, while the discretization where 

k^ = (5k)l^ /i = l,2,3,4 Z M = 0,l,...,iV-l (47) 

will be referred to as periodic. We emphasize that when we speak of periodic or antiperiodic 
case, we also imply that the corresponding reduction c(n) — > c(n) has been performed in a 
manner that is consistent with the remnant symmetries of the Riemann sum used. More pre- 
cisely, for the periodic discretization the Riemann sums I(j,n)[N] are completely symmetric 
with respect to indices n M and the reduction is thus as described in previous paragraph. 
However, for the antiperiodic case the Riemann sum is only symmetric with respect to ex- 
change of spatial components ni, n 2 and n%, and one thus has to keep track of larger set 
of coefficients c(n). It is only with this prescription that the exact correspondence between 
the finite and infinite volume situations is realized in case of antiperiodic discretization of 
the integral. One implication of this is that, for antiperiodic case at finite N, the time-like 
components of vectors and tensors will not be exactly equal to space-like components. How- 
ever, the difference (and recovery of full hypercubic symmetries of the infinite lattice) has 
to take place in the N — > 00 limit. In what follows, we will show results for both types of 
discretization. Perfect agreement has always been found as expected, and we thus provide 
the list of reduced coefficients c(n) of the periodic case only, which is more economic. 



3.4.1 The Results 

The first group of terms in Eq. (142]) is not gauge invariant and so constants i^ are expected 
to be zero. We have numerically evaluated these constants using the result for correspond- 
ing coefficients c(n) obtained from Mathematica. Fig. [2] shows the typical convergence of 
the associated Riemann sums i^[N] both in antiperiodic and periodic discretizations. For 
antiperiodic case both the time-like and the space-like case is shown to see that they are 
different at finite N. For iV — > 00 all three sequences are expected to approach the common 
limit which is apparently zero in this case. The convergence is again exponential in N. 
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The second and the third group of terms contain combinations that appear in (//, v fixed) 

tr c F^ u F^ = -tr c [A„, A v ] [A^ A u ] = 2 tr c (A fJi A fJi A u A u - A^A^A,) (48) 

and thus, in light of Conjecture 1, we anticipate that (3^ u = —a^ v . Also, due to hypercubic 
symmetries we must have that a^ u = a and (3^ v = (3 are independent of // ^ v. Since our 
results were obtained using a symbolic algebra software, verifying these relations is not only 
a consistency check for the expected result, but also an internal check of our programs. 

The set of reduced coefficients specifying constants and [3^ are given in Tables 1 
and 2 respectively (j = 9/2). In Fig. [3] (top) we show the typical convergence of Riemann 
sums for in both periodic and antiperiodic discretizations of the integral. For periodic 
case only ol^i is shown since full hypercubic symmetry in this case guarantees that, even 
at finite N, a^ u [N] is independent of //, v. Our results, obtained using symbolic algebra 
software, indeed comply with this. In the antiperiodic case we show both the representative 
of space-space combination (a^) and the representative of space-time combination (a^). 
All other possibilities for space-space and space-time combinations are exactly the same as 
the ones shown (even at finite N) as expected from the restricted hypercubic symmetry. The 
convergence of the three Riemann sums is exponential, as was in all the previous cases, and 
the insert shows that they approach a common non-zero limit whose value a can be easily 
extracted. The analogous results in case of (3^ v are shown in the bottom plot of Fig. [31 and 
one can immediately see that, indeed, (3 = —a. 

Combining all the results for the scalar spinorial component of D Q together, we have 
from Eq. (jUJ) that 

ti sc D Qj0 (U(a)) = tr sc A),o(I) + a 4 - J^tr + C(a 6 ) (49) 

implying that Conjecture 1 (see Eq.©) is indeed true for the class of constant fields. The 
associated proportionality constant is given by 

c S(p ; r ) = ^Efl _^ C S ( i \ = 0.01672926781 (50) 
4 V 19 / 

where the quoted value at p = 26/19 and r = 1 (only valid digits are given) was extracted 
from the data shown in Fig. [31 



4 Finite Volume and Non— Constant Fields 

As already mentioned in the previous section, we have analyzed the case of constant fields 
in the infinite volume in a way that can be directly mapped on to the situation in finite 
volume with specific boundary conditions. We now discuss this correspondence explicitly. 
Consider a 4-d hypercube in the continuum with side L p (in physical units). Superimposing 
a hypercubic lattice with iV 4 sites implies the classical lattice spacing a = such that 

L p = Na N (51) 
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Figure 3: Riemann sums for constants a^ u (top) and (3^ (bottom) of the scalar term (see 
Eq. fj42|) ). These results were computed at r = 1 and n — 0.19 (p = 26/19). 
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and the classical continuum limit is achieved as iV — > oo (ajv — > 0). Proper specification of 
the operator now involves also a choice of boundary conditions, and we will consider (see 
below) two standard cases that admit a diagonal Fourier-space representation for constant 
fields. In particular, we will refer to the situation with periodic boundary conditions in all 
directions as "periodic be", and to the case with periodic boundary conditions in spatial 
(/i = 1,2,3) directions and antiperiodic in time (/x = 4) as "antiperiodic be". The formulas 
([TO]) relating the direct and Fourier representations carry over to these finite volume cases 
except that the infinite sum is replaced by the finite sum, and the integral of the inverse 
transform is replaced by 



0, 



(27T)' 



J d A k e^-^ k O(k) — O n>m = el(n ~ m)k °^ ( 52 ) 



with sum only running over the associated N 4 discrete momenta. 

The relation to the infinite-volume case is based on the following considerations. 

(i) With periodic be, the iV 4 -dimensional Fourier space is spanned by plane waves with 
momenta given by Eq. ( 1471) . namely those involved in the periodic discretization (into N 4 
pieces) of the Brillouin zone discussed in case of infinite volume. Similarly, for antiperiodic 
be the associated momenta are given by Eq. (|40j) of the antiperiodic discretization. 

(ii) The form of Fourier representation D(k) of the overlap operator for the above finite 
volume setups (and constant fields) is identical to the infinite volume case but restricted to 
the corresponding discrete momenta. 

(Hi) Using (i), (ii) and comparing ( 1521 with the rule for evaluating Riemann sums in case 
of infinite volume (see e.g. Eq. (I39I) ) we can immediately see the following correspondence 

D$(a,N) = D${a)[N] (53) 

where the superscripts "inf" and "fin" refer to infinite and finite volume cases respectively. 
In other words, the matrix D^^a, N) (for finite volume setups discussed here) is equal to the 
corresponding Riemann sum (with N 4 terms) for the inverse Fourier transform representation 
of D™l(a). The same constant background is of course implicitly assumed in both cases. 

(iv) Riemann sums in question converge exponentially in N, as emphasized in previous 
sections, and we can conclude from (I5B1 that for large N 

D${a, N) = D$(a) + 0(eM~ C{a) N)) (54) 



0,0 v , 

with strictly positive C(a) (also at a = 0). 

From the last equation it follows that, in the case of constant background fields, the 
Taylor expansion (in a) of Dq o in infinite volume supplies an asymptotic expansion for the 
case of finite volume (with corections behaving as exp(— C(a)L p /a). Moreover, considering 
the continuum limit at fixed L p (so that = L p /N) we obtain that 

tr sc (D$(a N , N) - Z$(0, N)) tr Jtfgfta*) - D$(0)) 

lim ^ j — — <- = lim ^ — — — '- = c b (55) 

n^oo a% tr c Ff^Ffjo, n^oo a% tr c F^F^ 
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where c is the constant that we computed in the infinite volume. The above equation 
expresses the fact that the scalar component of D Q has the same classical limit in both 
finite and infinite volumes, as required and expected on the basis of locality. Analogous 
result obviously holds also for the tensor component. 

4.1 Non-Constant Fields 

Apart from the fact that the finite-volume setups are actually those relevant for practical 
lattice QCD calculations, their advantage for the present discussion is that one can evaluate 
the associated classical limit numerically. In other words, one can obtain the left-hand 
side of Eq. (1551) via direct numerical evaluation of Dq ™ for arbitrary classical background. 
Performing a sequence of such calculations with increasing N then allows us to infer the 
iV — > oo (ajy — > 0) limit of the ratio. For the case of constant fields, such calculations indeed 
exactly reproduce our results obtained via the expansion in classical lattice spacing. The 
utility of this approach however mainly lies in the fact that it allows us to evaluate classical 
limits also for non-constant backgrounds which were not included in our considerations up 
to this point. 

To perform such a calculation, we set L p — 1 and use arbitrarily selected classical back- 
grounds. One of the cases that we studied is specified by 

( sin27rt cos27rx/i\ / cos27T2:\ 
A^t, x,y,z)= 11 + — 1 a M + I 1 + ^ I a 3 _ M (56) 

where we use \i G {0, 1, 2, 3} and the coordinate-labeling correspondence <-> t, 1 x, 
2 <-» y, 3 <-> z for convenience. The constant field a M is specified in Appendix O 

For given N, we evaluated the ratio c s (N = 1/ajv) at 10 randomly selected points on 
the unit torus. The same points, specified in Appendix [Cj were used for each N. In Fig. @] 
(top) we plot the mean value of c s (N) over this sample with "error bars" representing the 
square root of the associated variance. The horizontal line shows the value of c s obtained 
via expansion in classical lattice spacing in infinite volume (1501) . As can be seen quite clearly, 
the mean approaches this predicted value for large N with variance shrinking toward zero 
at the same time. To guide the eye, we also included afittoiV = l/a7v dependence of c s in 
the form 

c S (N) = E c Hiv) (57) 

Note that this form is motivated by the fact that the odd powers of lattice spacing are not 
present in the expansion of L>o,o i n infinite volume. The above fitting form is neglecting the 
terms exponentially small in N, which are present, but their relative contribution decays 
very fast with increasing N. The standard fit shown in Fig. H] uses all points with iV > 12 
(treating the associated rooted variances as error bars) and K — 6. The coefficient Cq of 
the fit represents the estimate of c 5 , and the obtained value agrees with that of Eq. ([501 to 
relative precision of about two parts in 10 4 . 

To further check the robustness of the agreement found above (universality of the classical 
limit for constant and non-constant configurations), let us now examine a different procedure 
for estimating c s on background (|56l) . For any given point on the torus, consider evaluating 
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Figure 4: Evaluation of classical limit on a finite torus (L p = 1) for non-constant background 
field fl56|) . The meaning of the data is explained in the text. The overlap operator used is 
specified by r = 1 and k = 0.19 (p = 26/19). 
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the ratio (I55I) over some range of N (such as one shown in Fig. H] (top)), and then fitting 
this dependence to the form (!57|) with progressively increasing value of For each K, 

the coefficient Co = co(K) of the fit represents the estimate of c s . Since the expansion (JBTj) 
is asymptotically valid at large values of N, the accuracy of such estimates will depend on 
the fitting range chosen. In order to be very conservative, we use the whole range displayed 
in Fig. H] (top), i.e. all the points with 8 < N < 44. In Fig. H] (bottom) we show the average 
of Cq(K) so obtained over the same sample of random points that we used above. The 
"errorbars" on this plot are the square roots of variance over this sample, and are decaying 
rapidly with increasing K. This is expected since the classical limit has to be reproduced 
at each point individually and the estimates are expected to get better with increasing K. 
The horizontal line on the plot represents the value of c s quoted in Eq. (|50|) . As one can 
see quite clearly, c (K) settles approximately at this expected value. In fact, for the largest 
value of K (K = 9) the relative difference between c s and Co is about three parts in 10 5 . 

5 Values of Related Constants 

Given the results presented in Sees. El @] there is little doubt that the universal classical 
limit for the scalar spinorial part of -Do,o exists in the case of overlap operator, and that 
it is proportional to -^^(0)^^(0) pp. Indeed, we obtained fully consistent results showing 
this for finite and infinite volumes, for constant and non-constant fields, and for different 
boundary conditions in finite volume. Since matrix elements of D are local functions of the 
gauge field [5], the lattice operator 



represents a valid definition of tr cF^F^. For certain applications, such as those arising 
in the studies of QCD vacuum structure, it is necessary to deal with properly normalized 
operators in which case it is important that sufficiently precise values of c s and free-field 
subtraction constants are available. In this section, we provide some of this information both 
for the scalar and tensor cases. 

As is well known, the properties of the overlap Dirac operator, such as its range of 
locality [5], are quite sensitive to the value of the negative mass parameter p. Since various 
choices are currently being used in practice, we computed c s (p,r) and c T (p, r) for the range 
of p-values and the most commonly used case of r = 1. To obtain the results listed in 
Table El we used the expansion method described in Sec. El The last column lists the 
maximal discretization parameter N used in evaluating the associated Riemann sums. Note 
that one needs to use finer discretization (larger N) close to boundary values p = and 
p = 2 to achieve given precision. 

In Figure 0, we plot the p-dependence of both c s / p and c T / p. It is worth pointing out 
that in the tensor case the behavior appears to be almost exactly linear in p. However, 
closer inspection suggests that the deviations from exact linearity are in fact larger than the 

10 Note that all the points in the chosen range contribute with equal weight in such fit - there are no 
errorbars. Indeed, classical limit is not a "statistical" notion and has to be consistently reproduced at all 
non-singular points and for all classical configurations. 




(58) 
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Table 3: The proportionality constants c and c for various values of p at r = 1. 



estimated errors of calculated values. Nevertheless, the purely linear approximation is very 
precise in the range of p-values studied. 

Finally, we have computed the subtraction constants tr sc D ,o(I) for various sizes of a 
symmetric lattice with both periodic and antiperiodic boundary conditions. The results, 
summarized in Table HI were computed at p = 26/19 and r = 1. These are the values of 
the parameters used by the Kentucky group in the studies of QCD vacuum structure (see 
e.g. |3 H31 [21]) as well as in the studies of hadron spectroscopy (see e.g. [251 126]). In the 
space-time symmetric geometry considered here, the utility of these results is mainly for the 
former. To obtain the constants quoted in Table HI we directly evaluated tr sc D ,o(I) from 
the finite spectral sum over explicitly available Fourier modes. The calculation was done in 
double precision and we estimate that after rounding errors, these values are good to at least 
thirteen digits. Note that the convergence to the infinite volume value is very fast and that 
at iV = 40 the periodic and the antiperiodic cases already agree to 13 digits. 

6 Summary 

We have used analytic and numerical techniques to evaluate the classical limit of the gauge 
operator tr sc _D ,o(^) with D = D^'^ being the overlap Dirac matrix. As suggested 
on general grounds [lj, we found that the limit is proportional to tr c F^ V (Q)F^ V {Q) of the 

11 All the results presented in this manuscript were obtained in the Summer of 2006 and were partially 
discussed by the second author in Lattice 2006 talk that focused on coherent LQCD. In the Summer of 2007 
we were informed by David Adams [20] that he has just computed the constant c s (p, r) as well. His approach 
leads to an integral expression with integrand different from ours, but yielding identical final answers after 
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Figure 5: Dependence of rescaled constants c s (top) and c T (bottom) on the negative mass 
parameter p at r = 1. 
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N periodic be 



ant iper iodic be 



04 
08 
12 
16 
20 
24 
28 
32 
36 
40 



21.4358235894802 
21.3940949198520 
21.3931417096880 
21.3930879676987 
21.3930838480598 
21.3930834815890 
21.3930834459514 
21.3930834422695 
21.3930834418726 
21.3930834418356 



21.4010898463124 
21.3935062781772 
21.3931106079229 
21.3930856400191 
21.3930836424736 
21.3930834616044 
21.3930834438807 
21.3930834420471 
21.3930834418505 
21.3930834418306 



Table 4: The subtraction constants tr sc .Do,o(I) o n iV 4 lattice with periodic and antiperiodic 
(in time) boundary conditions. These results were calculated at p = 26/19 and r = 1. 

associated classical background after the subtraction of the free-field constant. Accordingly, 
the version of coherent LQCD where Tr D(U) serves as a basis for the gauge action [1] 
represents a valid regularization of QCD. In addition, this operator is expected to be useful 
as a natural partner to the overlap-based definition of pseudoscalar (topological) density both 
in studies of QCD vacuum structure and in standard applications of hadronic physics (such 
as calculations of glueball masses). In the former case it is of interest to explore the relation 
between the fundamental structure seen in topological density [121 131 [ZD [221 [231 121] , and 
the structure in scalar density that is expected to be visible if the overlap-based definition 
( 158]) is used. 

Techniques used here for the scalar case were applied in parallel also to the tensor com- 
ponent of Dq q- In our previous work [17] we have explicitly derived this classical limit 
(proportional to F^ u (0)) in a general setting, and the complete agreement found in this work 
serves as a valuable cross-check for both types of calculations. 

With possible practical applications in mind, we have computed the proportionality con- 
stants c and c T for a wide range of negative mass parameter p. One notable feature of the 
/^-dependence found is that c T (p) / p is surprisingly well described by a purely linear behavior 
in the region 0.2 < p < 1.8. 



integration. We hope that his derivation will be publicly available in due course. 



20 



A Alternative derivation to second order 



The expansion of the overlap operator for constant fields to the second order in the lattice 
spacing can be simplified as follows. We start by noticing that 

= (X^l B + S 2 , (59) 



where 



9X n T <9X dXj dX 



dk v dk v dk L 

Using the above formulas and the straightforward algebra one can show that 

1- v-^ . d „ 1/2 

'o > 



(60) 



lim — -£> = V^-I„B 

(61) 



lim = 

a^o aer p 

The derivative terms '^2 fl A fl -^f(k) vanish upon momentum integration and we have 

- p D 0fi (U(a)) = -A^I) + JJ^y, J d*k(- l -X»B-* /2 5^ + 0(a 3 ) (62) 
After expanding 5 2 and dropping the terms odd in k^ we get 

ifl^W) = i Do ,„(I) + £ (^E-- >< F ^ + ° (a3) (63) 

with £„„ defined in Eq, 
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B Reduced Coefficients 

The tables in this Appendix specify the reduced coefficients c(n) and powers j defining (via 
Eq. ( f46i) ) various constants appearing in the expansion of the scalar part of Z)o,o i n classical 
lattice spacing. As discussed in Sees. 13.41 and HJ these coefficients can also be used for 
evaluation of corresponding constants on a finite symmetric torus with periodic boundary 
conditions in all directions. 



m 


n 2 


™3 




c(n)/p 














48k 4 + 192r 2 K 4 + 1536k 6 - 1536r 2 K 6 - 10240r 4 K 6 + 12288k 8 - 73728r 2 K 8 + 122880r 4 K 8 


i 











-32rK 3 - 11527-k 5 - 2560r 3 K 5 - 18432rK 7 - 24576rV + 245760r 5 K 7 - 131072rK 9 + 










589824r 3 K 9 - 983040r 5 K 9 


i 


1 








576r 2 K 4 + 12416r 2 K 6 - 30720r 4 K 6 - 14336r 2 K 8 + 516096r 4 K 8 - 1474560r 6 K 8 


i 


1 


1 





-75264r 3 K 7 + 368640r 5 K 7 + 204800r 3 K 9 - 1179648r 5 K 9 + 1966080r 7 K 9 


i 


1 


1 


1 


-15360r 4 K 6 + 132096r 4 K 8 - 614400r 6 K 8 


2 











-192k 4 - 96r 2 K 4 - 5760k 6 + 15232r 2 K 6 + 20480r 4 K 6 - 43008k 8 + 321536r 2 K 8 - 










319488r 4 K 8 - 491520r 6 K 8 


2 


1 








64rK 3 + 4032rK 5 + 1216r 3 K 5 + 74496rK 7 - 232704r 3 K 7 + 290816rK 9 - 1503232r 3 K 9 + 










1179648r 5 K 9 + 2949120r 7 K 9 


2 


1 


1 





-1440r 2 K 4 - 13824r 2 K 6 - 42496r 4 K 6 - 198656r 2 K 8 + 952320r 4 K 8 - 1474560r 6 K 8 


2 


1 


1 


1 


6720r 3 K 5 - 18432r 3 K 7 + 242688r 5 K 7 + 77824r 3 K 9 - 774144r 5 K 9 + 1966080r 7 K 9 


2 


2 








240k 4 - 480r 2 K 4 + 8448k 6 - 25600r 2 K 6 - 23808r 4 K 6 + 77824k 8 - 561152r 2 K 8 + 861696r 4 K 8 + 










184320r 6 K 8 


2 


2 


1 





-5568rK 5 + 16896r 3 K 5 - 78336rK 7 + 170496r 3 K 7 + 344064r 5 K 7 - 282624rK 9 + 1835008r 3 K 9 - 










3336192r 5 K 9 + 737280r 7 K 9 


2 


2 


1 


1 


24576r 2 K 6 - 92928r 4 K 6 + 110592r 2 K 8 - 12288r 4 K 8 - 1081344r 6 K 8 


2 


2 


2 





-3840k 6 + 28416r 2 K 6 - 37888r 4 K 6 - 49152k 8 + 338944r 2 K 8 - 448512r 4 K 8 - 433152r 6 K 8 


2 


2 


2 


1 


23040rK 7 - 184320r 3 K 7 + 310272r 5 K 7 + 73728rK 9 - 518144r 3 K 9 + 516096r 5 K 9 + 1075200r 7 K 9 


2 


2 


2 


2 


7680k 8 - 67584r 2 K 8 + 241152r 4 K 8 - 224256r 6 K 8 


3 











576rK 5 + 64r 3 K 5 + 11520rK 7 + 25344r 3 K 7 - 122880r 5 K 7 + 184320rK 9 - 765952r 3 K 9 + 










786432r 5 K 9 + 327680r 7 K 9 


3 


1 








-288r 2 K 4 - 20224r 2 K 6 + 27136r 4 K 6 + 38912r 2 K 8 - 509952r 4 K 8 + 983040r 6 K 8 


3 


1 


1 





4800r 3 K 5 + 132096r 3 K 7 - 242688r 5 K 7 - 397312r 3 K 9 + 1306624r 5 K 9 - 983040r 7 K 9 


3 


1 


1 


1 


-15360r 4 K 6 - 153600r 4 K 8 + 307200r 6 K 8 


3 


2 








-2112rK 5 + 3776r 3 K 5 - 69888rK 7 + 116736r 3 K 7 + 76800r 5 K 7 - 540672rK 9 + 1982464r 3 K 9 - 










1628160r 5 K 9 - 1228800r 7 K 9 


3 


2 


1 





44544r 2 K 6 - 105472r 4 K 6 + 409600r 2 K 8 - 657408r 4 K 8 - 188416r 6 K 8 


3 


2 


1 


1 


-156672r 3 K 7 + 420864r 5 K 7 - 260096r 3 K 9 + 233472r 5 K 9 + 131072r 7 K 9 


3 


2 


2 





50688rK 7 - 238080r 3 K 7 + 258048r 5 K 7 + 432128rK 9 - 1789952r 3 K 9 + 1687552r 5 K 9 + 










792576r 7 K 9 


3 


2 


2 


1 


-211968r 2 K 8 + 1069056r 4 K 8 - 1422336r 6 K 8 


3 


2 


2 


2 


-73728rK 9 + 430080r 3 K 9 - 964608r 5 K 9 + 665600r 7 K 9 


3 


3 








8576r 2 K 6 - 11008r 4 K 6 - 18432r 2 K 8 + 132096r 4 K 8 - 184320r 6 K 8 


3 


3 


1 





-130560r 3 K 7 + 205824r 5 K 7 + 188416r 3 K 9 - 634880r 5 K 9 + 401408r 7 K 9 


3 


3 


1 


1 


297984r 4 K 8 - 528384r 6 K 8 


3 


3 


2 





-100352r 2 K 8 + 482304r 4 K 8 - 496640r 6 K 8 


3 


3 


2 


1 


237568r 3 K 9 - 1124352r 5 K 9 + 1382400r 7 K 9 


3 


3 


3 





-20480r 3 K 9 - 75776r 5 K 9 + 133120r 7 K 9 



Table 1 (part 1). Coefficients c(n) for constants a^ u of Eq. P2|) . In this case j = 9/2. 
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n 2 


n 3 




c{n)/p 


4 











-384k 6 - 3328r 2 K 6 + 3328r 4 K 6 - 10752k 8 - 28672r 2 K 8 - 109056r 4 K 8 + 266240r 6 K 8 


4 


1 








384rK 5 - 192r 3 K 5 - 3072rK 7 + 140544rV - 168960r 5 K 7 + 111616rK 9 - 188416r 3 K 9 + 
885760r 5 « 9 - 1351680r 7 K 9 


4 


1 


1 





-6912r 2 K 6 + 6912r 4 K 6 + 52224r 2 K 8 - 829440r 4 K 8 + 1019904r 6 K 8 


4 


1 


1 


1 


23040rV - 32256r 5 K 7 - 61440r 3 K 9 + 694272r 5 K 9 - 847872r 7 K 9 


4 


2 








2176r 2 K 6 - 3072r 4 K 6 - 9728k 8 + 251904r 2 K 8 - 430080r 4 K 8 + 64000r 6 K 8 


4 


2 


1 





-4608rK 7 - 39936r 3 K 7 + 70656r 5 K 7 - 26624rK 9 - 770048r 3 K 9 + 2027520r 5 K 9 - 997376r 7 K 9 


4 


2 


1 


1 


23040r 2 K 8 + 79872r 4 K 8 - 155136r 6 K 8 


4 


2 


2 





7680k 8 - 117760r 2 K 8 + 290304r 4 K 8 - 156672r 6 K 8 


4 


2 


2 


1 


3072rK 9 + 180224r 3 K 9 - 654336r 5 K 9 + 460800r 7 K 9 


4 


3 








19200rK 7 - 68352r 3 K 7 + 52224r 5 K 7 + 10240rK 9 - 33792r 3 K 9 - 182272r 5 K 9 + 332800r 7 K 9 


4 


3 


1 





-144384r 2 K 8 + 651264r 4 K 8 - 592896r 6 K 8 


4 


3 


1 


1 


125952r 3 K 9 - 712704r 5 K 9 + 721920r 7 K 9 


4 


3 


2 





-38912rK 9 + 417792r 3 K 9 - 950272r 5 K 9 + 563200r 7 K 9 


4 


4 








10240k 8 - 80384r 2 K 8 + 132096r 4 K 8 - 64512r 6 K 8 


4 


4 


1 





-34816rK 9 + 303104r 3 K 9 - 595968r 5 K 9 + 337920r 7 K 9 


5 











6912rK 7 - 768r 3 K 7 - 6144r 5 K 7 - 21504rK 9 + 131072r 3 K 9 + 7168r 5 K 9 - 188416r 7 K 9 


5 


1 








-3072r 2 K 6 + 2816r 4 K 6 - 25600r 2 K 8 - 181248r 4 K 8 + 237568r 6 K 8 


5 


1 


1 





41472r 3 K 7 - 41472r 5 K 7 - 63488r 3 K 9 + 706560r 5 K 9 - 765952r 7 K 9 


■5 


1 


1 


1 


-92160r 4 K 8 + 98304r 6 K 8 


■5 


2 








-9216rK 7 + 21504r 3 K 7 - 10752r 5 K 7 + 115712rK 9 - 458752r 3 K 9 + 489472r 5 K 9 - 64512r 7 K 9 


■5 


2 


1 





64512r 2 K 8 - 196608r 4 K 8 + 107520r 6 K 8 


■5 


2 


1 


1 


-64512r 3 K 9 + 251904r 5 K 9 - 162816r 7 K 9 


■5 


2 


2 





-52224rK 9 + 137216r 3 K 9 - 58368r 5 K 9 - 49152r 7 K 9 


5 


3 








-6656r 2 K 8 + 46080r 4 K 8 - 42496r 6 K 8 


5 


3 


1 





63488r 3 K 9 - 258048r 5 K 9 + 247808r 7 K 9 


5 


4 








-41984rK 9 + 156672r 3 K 9 - 189440r 5 K 9 + 76800r 7 K 9 


6 











4608k 8 - 25600r 2 K 8 + 24576r 4 K 8 - 3584r 6 K 8 


6 


1 








-2304rK 7 + 6912r 3 K 7 - 4608r 5 K 7 - 35840rK 9 + 112640r 3 K 9 - 10240r 5 K 9 - 66560r 7 K 9 


6 


1 


1 





19968r 2 K 8 - 64512r 4 K 8 + 44544r 6 K 8 


6 


1 


1 


1 


-21504r 3 K 9 + 73728r 5 K 9 - 52224r 7 K 9 


6 


2 








-7680k 8 + 42496r 2 K 8 - 58368r 4 K 8 + 23552r 6 K 8 


6 


2 


1 





34816rK 9 - 182272r 3 K 9 + 258048r 5 K 9 - 110592r 7 K 9 


c 




3 








ZJooZtk — 62464r k + 37888r k + 1024r k 


7 











-4096rK 9 + 17408r 3 K 9 - 22528r 5 K 9 + 9216r 7 K 9 


7 


1 








1536r 2 K 8 - 3072r 4 K 8 + 1536r 6 K 8 


7 


1 


1 





-9216r 3 K 9 + 18432r 5 K 9 - 9216r 7 K 9 


7 


2 








51207-K 9 - 23552r 3 K 9 + 31744r 5 K 9 - 13312r 7 K 9 


8 


1 








-1024rK 9 + 3072r 3 K 9 - 3072r 5 K 9 + 1024r 7 K 9 



Table 1 (part 2). Coefficients c{n) for constants a^ v of Eq. ( 1421) . In this case j = 9/2. 
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n 2 


n 3 


/7<4 


c\n)/ p 














-5120r 4 K 6 + 61440r 4 K 8 


l 











1024r 3 K 5 - 49152r 3 K 7 + 122880r 5 K 7 + 98304r 3 K 9 - 491520r 5 K 9 


l 


1 








4480r 2 K 6 - 24576r 4 K 6 - 71680r 2 K 8 + 589824r 4 K 8 - 737280r 6 K 8 


l 


1 


1 





-53760r 3 K 7 + 147456r 5 K 7 + 286720r 3 K 9 - 1179648r 5 K 9 + 983040r 7 K 9 


l 


1 


1 


1 


107520r 4 K 8 - 196608r 6 K 8 


2 











-96r 2 K 4 + 128r 2 K 6 + 7168r 4 K 6 - 45056r 2 K 8 - 49152r 4 K 8 - 245760r 6 K 8 


2 


1 








1856r 3 K 5 + 15360rK 7 + 29952r 3 K 7 - 55296r 5 K 7 - 40960rK 9 + 311296r 3 K 9 - 294912r 5 K 9 + 










1474560rV 


2 


1 


1 





-6656r 4 K 6 - 87040r 2 K 8 - 141312r 4 K 8 - 73728r 6 K 8 


2 


1 


1 


1 


-3072r 5 K 7 + 51200r 3 K 9 + 36864r 5 K 9 + 245760r 7 K 9 


2 


2 








-48k 4 + 96r 2 K 4 - 1536k 6 + 2176r 2 K 6 - 16128rV - 3328k 8 + 35840r 2 K 8 + 59904r 4 K 8 + 










258048r 6 K 8 


2 


2 


1 





960rK 5 - 1920r 3 K 5 + 18432rK 7 - 50688r 3 K 7 + 144384r 5 K 7 + 33792rK 9 - 28672r 3 K 9 + 










18432r 5 K 9 - 737280r 7 K 9 


2 


2 


1 


1 


-3840r 2 K 6 + 7680r 4 K 6 - 36864r 2 K 8 + 150528r 4 K 8 - 190464r 6 K 8 


2 


2 


2 





768k 6 - 5376r 2 K 6 + 10240r 4 K 6 + 12288k 8 - 80896r 2 K 8 + 43008r 4 K 8 - 248832r 6 K 8 


2 


2 


2 


1 


-4608rK 7 + 32256r 3 K 7 - 61440r 5 K 7 - 24576rK 9 + 149504r 3 K 9 - 184320r 5 K 9 + 399360r 7 K 9 


2 


2 


2 


2 


-1536k 8 + 12288r 2 K 8 - 33792r 4 K 8 + 7680O 6 k 8 


3 











-64r 3 K 5 + 60672r 3 K 7 - 67584r 5 K 7 - 57344r 3 K 9 + 294912r 5 K 9 + 163840r 7 K 9 


3 


1 








-8960r 2 K 6 + 9728r 4 K 6 + 117760r 2 K 8 - 964608r 4 K 8 + 712704r 6 K 8 


3 


1 


1 





107520r 3 K 7 - 107520r 5 K 7 - 542720r 3 K 9 + 2330624r 5 K 9 - 1228800r 7 K 9 


3 


1 


1 


1 


-215040r 4 K 8 + 208896r 6 K 8 


3 


2 








960rK 5 - 896r 3 K 5 + 3072rK 7 - 70656r 3 K 7 + 98304r 5 K 7 + 64512rK 9 - 231424r 3 K 9 + 










374784r 5 K 9 - 983040r 7 K 9 


3 


2 


1 





-15360r 2 K 6 + 14336r 4 K 6 - 96256r 2 K 8 + 546816r 4 K 8 - 524288r 6 K 8 


3 


2 


1 


1 


46080r 3 K 7 - 43008r 5 K 7 + 167936r 3 K 9 - 528384r 5 K 9 + 262144r 7 K 9 


3 


2 


2 





-13824rK 7 + 105984r 3 K 7 - 104448r 5 K 7 - 94208rK 9 + 296960r 3 K 9 - 446464r 5 K 9 + 










878592r 7 K 9 


3 


2 


2 


1 


55296r 2 K 8 - 423936r 4 K 8 + 417792r 6 K 8 


3 


2 


2 


2 


12288rK 9 - 73728r 3 K 9 + 282624r 5 K 9 - 352256rV 


3 


n 

6 








p a a 2 6 cio 4 6 onn/1/1 2 8 i ooTnon 4 8 oaotci 6 8 

64Ur k — 51zr k — 82944?" k + 33792l)r k — 2U2752r k 


3 


3 


1 





-7680r 3 K 7 + 6144r 5 K 7 + 475136r 3 K 9 - 1404928r 5 K 9 + 679936r 7 K 9 


3 


3 


1 


1 


15360r 4 K 8 - 12288r 6 K 8 


3 


3 


2 





91136r 2 K 8 - 353280r 4 K 8 + 229376r 6 K 8 


3 


3 


2 


1 


-182272r 3 K 9 + 706560r 5 K 9 - 458752r 7 K 9 


3 


3 


3 





-108544r 3 K 9 + 260096r 5 K 9 - 114688r 7 K 9 



Table 2 (part 1). Coefficients c(n) for constants (3^ u of Eq. fT42l . In this case j = 9/2. 
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n 2 


n 3 


// 1 


c[n)/ p 


4 











-2432r 2 K 6 + 1792r 4 K 6 + 37888r 2 K 8 - 167424r 4 K 8 + 151552r 6 K 8 


4 


1 








-153607-K 7 + 85248r 3 K 7 - 56832r 5 K 7 + 51200rK 9 - 456704r 3 K 9 + 1276928r 5 K 9 - 
860160r 7 K 9 


4 


1 


1 





122880r 2 K 8 - 506880r 4 K 8 + 325632r 6 K 8 


4 


1 


1 


1 


-122880r 3 K 9 + 448512r 5 K 9 - 282624r 7 K 9 


4 


2 








768k 6 - 2176r 2 K 6 + 2048r 4 K 6 - 5632k 8 + 1536r 2 K 8 + 61440r 4 K 8 - 134656r 6 K 8 


4 


2 


1 





-9216rK 7 + 26112r 3 K 7 - 24576r 5 K 7 + 2048rK 9 - 22528r 3 K 9 - 110592r 5 K 9 + 309248r 7 K 9 


4 


2 


1 


1 


18432r 2 K 8 - 52224r 4 K 8 + 49152r 6 K 8 


4 


2 


2 





-7680k 8 + 71680r 2 K 8 - 175104r 4 K 8 + 156672r 6 K 8 


4 


2 


2 


1 


15360rK 9 - 143360r 3 K 9 + 350208r 5 K 9 - 313344r 7 K 9 


4 


3 








10752rK 7 - 28416rV + 12288r 5 K 7 - 50176rK 9 + 347136r 3 K 9 - 727040r 5 K 9 + 404480rV 


4 


3 


1 





-86016r 2 K 8 + 227328r 4 K 8 - 98304r 6 K 8 


4 


3 


1 


1 


86016r 3 K 9 - 227328r 5 K 9 + 98304r 7 K 9 


4 


3 


2 





51200rK 9 - 245760r 3 K 9 + 495616r 5 K 9 - 331776r 7 K 9 


4 


4 








7424k 8 - 44032r 2 K 8 + 62976r 4 K 8 - 21504r 6 K 8 


4 


4 


1 





-29696rK 9 + 176128r 3 K 9 - 251904r 5 K 9 + 86016r 7 K 9 


5 











2304r 3 K 7 - 1536r 5 K 7 - 20480r 3 K 9 + 91136r 5 K 9 - 106496r 7 K 9 


5 


1 








25600r 2 K 8 - 86016r 4 K 8 + 51200r 6 K 8 


5 


1 


1 





-102400r 3 K 9 + 288768r 5 K 9 - 167936r 7 K 9 


5 


2 








-4608rK 7 + 9984r 3 K 7 - 6144r 5 K 7 + 1024rK 9 - 11264r 3 K 9 - 22528r 5 K 9 + 113664r 7 K 9 


5 


2 


1 





36864r 2 K 8 - 79872r 4 K 8 + 49152r 6 K 8 


5 


2 


1 


1 


-36864r 3 K 9 + 79872r 5 K 9 - 49152r 7 K 9 


5 


2 


2 





15360rK 9 - 81920r 3 K 9 + 141312r 5 K 9 - 116736r 7 K 9 


5 


3 








-7168r 2 K 8 + 18432r 4 K 8 - 8192r 6 K 8 


5 


3 


1 





28672r 3 K 9 - 73728r 5 K 9 + 32768r 7 K 9 


5 


4 








-19456rK 9 + 67584r 3 K 9 - 71680r 5 K 9 + 20480r 7 K 9 


6 











-6656r 2 K 8 + 12288r 4 K 8 - 5632r 6 K 8 





1 


u 


("1 

u 


— l(J24l)rK + 44L)32r k — olzUUr k + 174l)or k 


6 


2 








-1536k 8 + 17408r 2 K 8 - 29184r 4 K 8 + 13312r 6 K 8 


6 


2 


1 





6144rK 9 - 69632r 3 K 9 + 116736r 5 K 9 - 53248r 7 K 9 


6 


3 








13312rK 9 - 38912r 3 K 9 + 35840r 5 K 9 - 10240r 7 K 9 


7 











7168r 3 K 9 - 14336r 5 K 9 + 7168r 7 K 9 


7 


2 








3072rK 9 - 16384r 3 K 9 + 23552r 5 K 9 - 10240r 7 K 9 



Table 2 (part 2). Coefficients c{n) for constants (3^ of Eq. ( 1421) . In this case j = 9/2. 
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C Sample Constant Field and Sample Points 

The constant field a M used in numerical calculations of Sec. 14.11 (see Eq. fl56l) ) is specified by 



0.333027 
a = ( 0.438189 - 0.318128i 
0.374405 - 0.222978?' 



0.438189 + 0.318128? 

-0.154498 
0.378804 - 0.22751? 



0.374405 + 0.222978?' 
0.378804 + 0.2275 li 
-0.178529 



0.230383 0.202411 + 0.030092?' 0.384110 - 0.364239?' ' 

a = | 0.202411 - 0.030092?- -0.162558 0.448044- 0.138645?- 

0.384110 + 0.364239? 0.448044 + 0.138645? -0.067824 



-0.254314 
a 2 = \ 0.490871 + 0.296415z 
0.709724 + 0.352729i 



0.490871 - 0.296415? 

0.330503 
0.970160 - 0.109804i 



0.709724 - 0.352729i 
0.970160 + 0.109804?- 
-0.076189 



-0.055155 
(A 3 = \ 0.502826 - 0.023544?- 
0.435721 + 0.074859? 



0.502826 + 0.023544?- 

0.159857 
0.581485 - 0.032500?- 



0.435721 - 0.074859?- 
0.581485 + 0.032500? 
-0.104703 



The classical limit on non-constant backgrounds was evaluated at 
generated points on the unit symmetric torus. 



x 

0.103172983143503 
0.248194925172178 
0.174820791004537 
0.450532157114747 
0.416208702090845 
0.295893446617454 
0.192720463473952 
0.944525538301774 
0.769704747297237 
0.319172590182490 



y 

0.9458719000957324 
0.0893431271786154 
0.7529966955146729 
0.7958101168552333 
0.6286435317521117 
0.0590763301619042 
0.1132044300661718 
0.0238613028875564 
0.2708646073728835 
0.4750544905176503 



0.956116383180309 
0.336764036419274 
0.984186866641149 
0.274186255621803 
0.459165616384016 
0.678511201131179 
0.722394817950870 
0.385630781531596 
0.401443914890447 
0.127257659268644 



the following randomly 



t 

0.438505993255755 
0.118409829321844 
0.265965702089591 
0.321281009101127 
0.356014499573743 
0.781239724758651 
0.342733731502896 
0.224323902181052 
0.958358200091461 
0.637077190990334 
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